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The  double  sheaths  formed  by  plasmas  that  contain  thermal  electrons , as  well  as  positive  and 
negative  ions  has  been  investigated . The  negative  ion  beam  that  can  be  extracted  under  these  conditions 
is  given  as  a function  of  the  electric  potential , the  density  ratio  and  the  Langmuir  limit . The  result  is 
useful  when  it  comes  to  extract  ion  beams  from  negative  ion  sources . 


1.  Introduction 

Double  sheaths  set  up  by  electron  emission  [1] 
and  a counter-stream  [2]  can  be  tailored  to  give 
optimum  extraction  of  beams  with  regard  to  the 
Langmuir  limit.  Double  sheaths  containing 
negative  ions  are  scarce  [3].  In  a previous  work 
the  effect  of  negative  ions  within  a plasma  upon 
electron  beam  extraction  has  been  investigated 
[4].  This  paper  deals  with  double  sheaths  whereby 
the  plasmas  on  both  sides  do  contain  negative 
ions.  This  is  of  practical  interest  if  one  wants  to 
optimize  negative  ion  beams  from  an  NBI 
source. 

2.  Model 

Consider  a double  layer  formed  by  two 
plasmas,  both  containing  negative  ions.  From  the 
lower  potential  side  (potential  F=0),  thermal 
electrons  and  negative  ions  flow  with  current 
density  Jc  and  Jn.  Particle  densities  nc](V),  nn \(V) 
are  given  by 

»c.=  7 »„■  , (1) 

ef2e(V  +V€)lm  e^2e(V  + Vn)lMn 

where  Vc  and  Vn  are  the  potentials  equivalent  to 
the  initial  energy,  m>  Mp,  Mn,  are  electron,  positive 
and  negative  ion  masses.  Positive  ions  get 
reflected  and  their  density  np \(V)  within  the  sheath 
is 

npl(V)  = Np\txp(-eV/KTp),  (2) 

where  Np\  is  the  boundary  value.  From  the  high 
plasma  potential  side  ( V~ Va)  positive  ions  flow 
with  current  density  Jp,  such  that  the  density 
becomes 

np2(V)  = JP /epe(Va-V+V0)/Mp  , (3) 

whereby  the  initial  velocity  is  given  by  vG 
according  eV0=  Mpvfl2.  Electrons  and  negative 
ions  are  reflected.  Their  densities  are 

ne2  = Ne2  exp [-e(Va  - V)lKTe] , 
nn2  = N/i2 exp[-e(K0  - V)l kT}J]  , (4) 


where  Tc , Tn,  Tp , are  the  temperatures  and  N& , N P2, 
/Vn2  are  the  densities  at  the  boundary. 

3.  Formulation 

The  right  hand  side  of  Poisson’s  equation  is 
LI  iff  SQ  ~~  ^//l  w p\  ~ ne  2 — nn2  + n ^2  > (5) 

where  £0  is  the  dielectric  constant  in  vacuum. 

Introducing  the  following  normalization, 

£=x/Ad  ,XD=(e0KT>JNple2)'r\  ri=eV/KTc,  t]  = 
eVjKTc,  i)n=eVjKT',  vcl=Nc[/Nph  vnl=NJNpi, 
Ve2=Nc2/NP2,  vn2=jVn2/A/p2,  q=Np2/ Np\,  y-a=TJTm 
y p TJ  Tp,  L JJJa,  jn  qi  jp-JpOp/Jo,  where 
OH MJm)  , <V=  (Mp/m)'12,  J0=Wple(2Krc/m)l/2, 
Poisson's  equation  now  reads 

d2rj/d<f =-ep/£0 . P=N/Npl  (6) 

P = — — + exp(-y  77)  + - =JjL 

yjV+rie  yjTl  + TJ,,  V?7a-P  + P0 

-qve2  exp[-(rja  -ri)]-qv„2  exp[-y„  (rja  -77)] , (7) 
and  the  conditions  at  the  plasma  boundaries  are 

1)  Quasi  neutrality  : p=  0 at  77=0,  7]a. 

2)  Zero  derivative  of  space  charge  density 

dp/d£=0  at  77=0,  rja. 

3)  Zero  electric  field:  dp/dl^O  at  77=0,  rja. 
Multiplying  d77/d£  on  both  sides  of  (6)  and 

integrating,  we  obtain  the  stress  (d?7/d£)2/2,  which 
automatically  satisfies  the  condition  drj/d^O  at 
77=  77a.  At  77=0,  then  we  obtain 

0 = 2je(%U2 --Jva  +%)+V„{n„U2 -Ll +n„ ) 

+ [1  - e\p(-yp7]a )]/yp  + 2jp (^77,, +77,,  - ) 

+qve2{exp{-ria)-l}  + qv„2{exp(-ynria)-l}/yn . (8) 
The  double  layer  does  form,  if  eqn.(8) 
together  with  the  boundary  conditions  is  satisfied. 

4.  Double  sheath  with  high  potential 


Putting  exp(-77a)«l,  exp(-yp7ja)«l,  exp(-ye77a) 
«1  and  exp(-yn77a)«l  and  rjc=rjni  we  get 


Jc  Jn  _ j J P 

„ 1/2  LTTZ  ’ 


(9) 
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Jc  + Jn  _ ^ 

Jla+le  yl% 


-qfra  +v,»2). 


Jc  + ■/« 
3/2 

2'?,. 


-r»+- 


2(a,+fU 


3/2 


= 0, 


(10) 

(N) 


2(na+nc)1'2 


jn 

3/2 

2/)o 


= 9(^,2 + r„v„2), 


(12) 


i+y‘p  /('L+'Jo  )l/2 
2 Yp-jpKVa+Tlo) 


(13) 


2*7 (^t'2  y«^/i2  ) *“  (7  c*  7/i  ) ^ a 


(14) 


pn  and  p0  given  by  (13)  and  (14)  correspond  to 
the  Bohm  criterion  determining  the  ion  fluxes.  pe 
(or  pn)  increases  from  1/(2 yp)  as  pa  is  decreased  or 
yp  is  increased.  pc  increases  from  1/(2 yn)  for  the 
negative  ion-rich  plasma  as  pP  is  decreased  or  jc 
and  yn  are  increased.  p0  tends  to  1/2  of  the  simple 
sheath  if yc=yn-0  and  77a>>  1 . 

From  (8),  we  obtain  the  Langmuir  limit  as 

Ue  + Jn  )(A, +'?„  ~§h)~  jp(<Jria+n„  ~ VT ) 


= l/y/,-//(v(.2+v„2/y„).  (15) 

If  17,  is  extremely  large,  the  limit  is  modified  as 

J c Jn  = j p • (1^) 

Without  yn,  we  recover  Jc=(Mp/m)u2Jp.  Without 
yp,  we  obtain  Jn=(AYp/A7n)1/2yp,  i.e.  the  negative  ion 
beam  is  governed  by  the  positive  ion  beam,  and 
determined  by  the  square  root  of  the  mass  ratio. 


5.  Calculations  and  Results 


We  have  five  equations  to  determine  the 
unknowns  p0,  pc , yc+  yn=ycn,  yp  and  q from  the 
known  values  pa,  vc2,  vn2,  yn  and  yp.  We  define 
jcn~ yc+yn-  This  nonlinear  system  has  been  solved 

by 


Fig.  1.  7c n andycn/yp  vs  for  pa  =10  and  vc2=0.3. 


iteration  starting  with  p0=  1/2,  pn=I/(2y)  to  get  ycn, 
7P  and  q.  Afterwards  new  values  for  p0  and  pn  , 
e.t.c  were  determined. 

In  Fig.l  for  pa=l 0,  yn=yp=10,  the  beam  density 
y‘cn  has  a maximum  at  a vn2  while  pa  is  small  but 
becomes  almost  constant  for  higher  pa.  This 
behavior  holds  in  a wide  range. 

Fig.  2.  p,  7 f and  p vs  £ for  p3=  1 0,  yn=yp=  1 0. 


In  Fig.  2 showing  the  profile  of  the  double  sheath, 
p is  not  symmetric  as  in  the  beam-free 
case,  showing  a narrow  dip  near  the  low  potential 
side  and  a small  bump  on  the  higher  side,  though 
it  was  confirmed  that  the  integral  of  p over  the 
entire  region  is  0.  For  vn2=0.1,  p shows  a wavy 
structure  near  pa,  suggesting  a triple  layer.  As  vn2 
is  increased  above  0.1,  p,  p'  and  p resemble  those 
of  vn2=0.  Therefore,  except  the  cases  with  a triple 
layer,  the  sheath  thickness  6 (normalized  by  AD) 
can  be  defined  by  the  distance  between  0.01  pa 
and  0.99 pa.  8 deviates  from  p4/'Vy'cnl/2,  as  pa  is 
decreased. 

6.  Conclusion 

The  double  sheath  theory  has  been  extended  to 
the  case  when  a beam  of  negatively  charged  ions 
is  extracted  from  a plasma  containing  such 
particles.  Results  are  presented  for  the  current 
densities  of  positively  and  negatively  charged 
species,  profiles  of  space  charge,  electric  field  and 
potential  within  the  sheath.  The  potentials  at  the 
sheath  edge  will  help  to  better  understand  the 
extraction  of  negative  ions  from  negative  ion 
sources. 
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